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Abstract 

We investigate tunneling properties of collective excitations in the ferromagnetic phase of a spin-1 
spinor Bose-Einstein condensate (BEC). In addition to the Bogoliubov mode, this superfluid phase 
has two spin excitations, namely, the gapless transverse spin wave and the quadrupolar mode with 
a finite excitation gap. In the mean-field theory at T = 0, we examine how these collective modes 
tunnel through a barrier potential that couples to the local density of particles. In the presence of 
supercurrent with a finite momentum q, while the Bogoliubov mode shows the so-called anomalous 
tunneling behavior (which is characterized by perfect transmission) in the low energy limit, the 
transverse spin-wave transmits perfectly only when the momentum k of this mode coincides with 
ibg. At A; = itg, the wave function of this spin wave has the same form as the condensate wave 
function in the current carrying state, so that the mechanism of this perfect transmission is found 
to be the same as tunneling of supercurrent. Using this fact, the perfect transmission of the 
spin wave is proved for a generic barrier potential. We show that such perfect transmission does 
not occur in the quadrupolar mode. Further, we consider the effects of potentials breaking U(l) 
and spin rotation symmetries on the transmission properties of excitations. Our results would be 
useful for understanding excitation properties of spinor BECs, as well as the anomalous tunneling 
phenomenon in Bose superfluids. 

PACS numbers: 03.75.Lm,03.75.Mn,75.30.Ds,75.40.Gb 
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I. INTRODUCTION 



Recently, an anomalous tunneling phenomenon has been extensively discussed in the 
field of Bose gas superfiuids [l-16|. In this phenomenon, the Bogoliubov mode tunnels 
through a barrier without reflection in the low energy limit. This tunneling property is 
quite different from that of an ordinary single particle in quantum mechanics, where perfect 
reflection occurs in the low energy limit. The anomalous tunneling phenomenon occurs even 
in the presence of a finite superfiow, except in the critical supercurrent state. In the critical 
current state, tunneling of the Bogoliubov mode is accompanied by a finite refiection in the 



low energy limit js, 8,ll0|. Since the Bogoliubov mode is a collective Nambu-Goldstone (NG) 



mode associated with the broken U(l) gauge symmetry in superfiuid phases, the anomalous 
tunneling phenomenon may be considered as a fundamental property that Bose superfiuids 
generally have. 

While the Bogoliubov mode dominates over low energy properties of a (spinless) scalar 
Bose-Einstein condensate (BEG), a spinor BEG also has other collective modes associated 
with spin degrees of freedom. It is an interesting problem whether or not such spin wave 
modes also exhibit anomalous tunneling behaviors as in the case of the Bogoliubov mode. In 
particular, in the ferromagnetic spinor BEG 
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18| . in addition to the ordinary Bogoliubov 
mode, there is a gapless transverse spin wave mode with a quadratic dispersion E (x 
(where p is the momentum of the spin wave), which is similar to the magnon in a ferromag- 



net 



19|. Thus, using the ferromagnetic spinor BEG, one can conveniently study physical 



properties of both the type-I Nambu-Goldstone mode (Bogoliubov mode) characterized by 
E oc and type-II Nambu-Goldstone mode (transverse spin wave) characterized by 

E oc p^", where n is an integer 2q|. We also note that the ferromagnetic phase of a spin-1 
BEG has the other spin mode, called the quadrupolar mode, with a finite excitation gap. 
Thus, this superfiuid phase is a very useful system to examine how the detailed excitation 
spectra of collective excitations affect their tunneling properties. 

In this paper, we investigate tunneling properties of low energy collective excitations in 
the ferromagnetic phase of a spin-1 spinor BEG at T = 0. In a previous paper 
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16|, 



two of the authors showed the anomalous tunneling behavior of the 



in the absence of superfiow. In this paper, we extend Refs. 



12 



14- 



;ransverse spin wave 



16| to the case with 



a finite superfiow. Within the framework of the mean-field theory for the spin-1 BEG at 
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T = 0, we determine the spatial variation of condensate wave functions around the barrier 
by solving the Gross-Pitaevskii (GP) equation. We then clarify tunneling properties of 
collective excitations by solving the Bogoliubov equations. We show that the low-energy 
tunneling properties of spin wave excitations are very different from those of the Bogoliubov 
mode in the presence of a finite supercurrent. In particular, perfect transmission of the 
transverse spin wave is shown to occur through a generic potential barrier, not in the low 
momentum limit, but in the case when the magnitude of spin wave momentum \k\ and that 
of the momentum of the supercurrent |g| are equal. (Note that perfect transmission of the 
Bogoliubov mode always occurs in the zero momentum limit {k — )■ 0), irrespective of the 
value of g.) 

This paper is organized as follows. In Sec. [TTl we explain our formulation for tunneling of 
collective excitations through a barrier in current carrying states of a spin-1 ferromagnetic 
spinor BEG. We discuss tunneling properties of the transverse spin wave and the quadrupolar 
mode in Sees. IIIII and IIVI respectively. In Sec. |Vl we examine how the symmetry of a 
barrier potential affects the anomalous tunneling properties of the Bogoliubov mode and 
the transverse spin wave. 



II. MEAN-FIELD THEORY OF SPIN-1 BEC 



We consider the supercurrent state of a ferromagnetic spin-1 spinor BEG with a barrier 
potential put around x = 0, as schematically shown in Fig. [H Since the y- and 2;-direction 
are not important in this paper, we simply assume a one-dimensional system. The spin 
quantization axis is taken parallel to the 2;-direction, and we introduce a three-component 
condensate wave function 



(1) 



where is the component of magnetic sublevel j = Sz (±1 and 0). As usual, the con- 

densate wave function ^{x,t) obeys the GP equation obtained from the variational principle 
of the action 



dt / dxC. 



(2) 
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FIG. 1: (color online) Model of a spin-1 ferromagnetic spinor BEC with a potential barrier V{x) 
put around x = 0. The ferromagnetic condensate wave function is described hj Sz = 1 spin 
component, and a finite supercurrent flows in the j;-direction. In this paper, a collective mode 
is injected from x = — oo. For simplicity, we assume a one-dimensional system and ignore the 
unimportant y- and z-direction. 



Here, £ = ih^\x,t)d^ 



density T-L is given by 



X, t) — ^{x^t) is the Lagrangian density, where the Hamiltonian 
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18| 



^^^\x,t)dl^{x,t) + V{x)p{x,t) + yp'(x,t) + ^F2(x,t) - (7PBSF,(x,t). (3) 



^1 1:^2/ 



Here, V{x) is a barrier potential put around x = in Fig. [H Cq = 47r^^(ao + 2a2)/(3m) and 
Ci = ATih^{a2 — ao)/(3r?7,) describe, respectively, coupling constants of a spin-independent 
and spin- dependent interactions 18|, where ag and 02 are the s-wave scattering lengths be- 
tween two atoms in the total spin S* = state and S = 2 state, respectively. p(x, t) = 
{x,t)^{x,t) is the particle density, and F = ^{F^, Fy, F^) is the spin density, where 
Fj{x,t) = {x,t)Sj^{x,t). Since we are taking the spin quantization axis parallel to the 
z-direction, the S = 1 spin matrices Sj{j = x, y, z) are given by 
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The last term in Eq. ([3]) describes effects of a magnetic field B applied parallel to the 2;-axis, 
where g and /ib are the Lande's g factor and Bohr magneton, respectively. We take B to be 
non-negative without loss of generality. 

Taking the variation of the action / in Eq. ([2]) with respect to ^'^{x,t), we obtain the 
time- dependent GP equation, 

a$(x,t) _ 



ih- 



dt 



\ 







V2 



F+ h{x, t) - ciF^ + gfisB 



^x,t). (5) 



Here, h{x,t) = —h'^d'^/{2m) + V{x) + cop{x,t), and F± = i: iFy. The time-independent 
GP equation for stationary states is then obtained by setting = e~*^J*/''$j(x) under 



the condition 2/io 
/ 



/ii + [21]. In this paper, taking /ij-i 

Cl 



h{x) + CiF^ - guBB -^F. 

Cl 

V2 



h{x) 





V2 



\ 



fx 



21 



22 1, we obtain 



v 







—;=F+ h{x) - ciF^ + gfjiBB 
v2 / 



^ix) = 0, 



(6) 



where h{x) = —f^d"^/ {2m) — fi + V{x) + Cop(x). 

The ferromagnetic state is always realized as the ground state when the spin-spin inter- 
action is ferromagnetic (ci < 0) and B = 0. This region in ci-B plane becomes wider as 
Cl < gii^B/ Pq. In the ferromagnetic state, all the Bose atoms occupy the Sz = l state as 



'^$+i(x)^ 



<^{x) 



\ 








(7) 



In the uniform system {y{x) = 0), the Sz = 1 component in the supercurrent state has the 
form <l>_|_i(x) = y/poe^'^^, and the chemical potential p is given hj p = c+po + h^q^/{2m) — 
gpsB. Here, po is the particle density of the uniform Bose gas, and c+ = Cq + Ci is the 
coupling constant of the interaction between the atoms in the 8^ = 1 state. In addition, as 
will be shown soon later, the other spin components, $0(2^) and are not induced by 

the barrier V{x). Using these results, one finds that the GP equation in the ferromagnetic 
state reduces to 



+ V[x) - c+po - h c+p{x) 



2m dx"^ 



2m 



$+i(x) = 0. 



(8) 
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FIG. 2: (color online) Calculated condensate wave function ^j^i{x) = |$_|_i(x)|e*'''+i*^^-' in the 
supercurrent state of a ferromagnetic spin-1 spinor BEC at T = for the barrier potential in 
Eq. ^ with Vq = 20+ po- Since ^o{x) and <I>_i(a;) are absent everywhere, we do not show them 
here, (a) amphtude |<I>_|_i(x)|. (b) phase 99+1 (x). In panel (b), we set 93+1(0) = 0. The inset shows 
the momentum q of the supercurrent as a function of the phase difference A99. The red point 
{q = 0.0574,^"^ = Qc) in the inset shows the critical current state. 

We note that Eq. ([8]) has the same form as the GP equation for the scalar BEC, when we 
regard c+ as the coupling constant of an interaction between spinless Bose atoms. 

Figure [2] shows the spatial variation of the condensate wave function $+i(a;) = 
|$+i(x)|e*^+i(^). In this fi gure, we consider the potential barrier 

V{x) = Koe-^'/«^ (9) 

where ^ = h/ ^mc^p^ is the healing length. Although we have numerically solved the 
GP equation ([6]) to include the possibility of finite $0(2;) and $_i(s) components around 
X = 0, these components have turned out to be absent. Thus, one may safely use the 
simpler GP equation ([8]). As expected. Fig. [2] (a) shows that the condensate wave function 
is suppressed around the potential barrier. In addition, panel (b) shows that the phase 
99+1 (x) of the condensate wave function spatially varies. This behavior reflects the presence 
of a finite superfiow. Far away from the barrier, the phase factor reduces to the form 
V3+i(a; — > ±00) = gx + C± The phase difference = C+ — C_ can be conveniently 
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calculated from 



23| 



Ay? = g 



dx 



(10) 



The inset of Fig. [2] (a) shows the momentum q of the supercurrent as a function of A(f. As 
in the case of the scalar BEC, the upper critical momentum exists at Aip ~ 7r/2 (the red 
point in the inset). 

Once the condensate wave function $(x) is determined, collective excitations are con- 
veniently obtained by considering small fluctuations around ^{x). Setting <l>(x,t) = 
<l>(x) + 0(x, t) (0 = *(0+i, 00) 0-i))) cind substituting it into Eq. ([5]), one obtains the Bogoli- 
ubov equations [2^ 



ih- 



')±1 



dt 



h{x) T QPbB + i?it±i + Ci(±F. + 1$ 



+ p±i0o + i?Si0s + <t;±i0ti + (^Vi,-i + ci'^>o)0;i, 

= (h{x) + cip(x) + i?Jo^) 00 + (co$o + 2ci<l>+i$_i)0S 

+ $:(i^*0, + <.V*), 

i=±l 



(11) 



(12) 



where P±i = (co + Ci)$*<l>±i + 2ci$;i<l>o, = (co±ci)$,$,, and Pjf = (cq ±Ci)<l>,*$,. In 
obtaining Eqs. ( fTTj) and ( |T2i) . we have retained terms up to O(0(x)). In the ferromagnetic 
case, Eqs. (|TT1) and (|T2l) reduce to 



^+1 



at 



ih- 



dt 



ih- 



dt 



2m dx'^ 

h^d^ 

2m dx"^ 

d^ 



2J1 



^V[x) 
\V{x) 



c+Po 
c+po 



hrq 
2m 

2m 
h^q^ 



+ 2c+p(x) 



0+i + c+$^i(x)0;i, 



+ c+p(x) + g^isB 



+ V{x) - c+po - ir^ + (c+ - 2ci)p(x) + 2^pb5 



(13) 
(14) 

(15) 



2m dx"^ ' ^ ' " 2m 
Since Eqs. f|T3|) -f fT5|) are independent, we can treat each component 0j(x) separately. Equa- 
tion (fT3l) has the same form as the ordinary (time-dependent) Bogoliubov equation in the 
scalar BEC [24j. We thus see that the mode described by 0+i(a;) is the Bogoliubov mode in 
the present ferromagnetic spinor BEC and exhibits the anomalous tunneling phenomenon 
as in the case of the scalar BEC jl-^. Namely, the Bogoliubov mode in the low energy limit 
always tunnels through the barrier without reflection in the supercurrent state 0, 8], except 
in the critical current state. In the critical supercurrent state, the Bogoliubov mode tunnels 



through the barrier with a finite reflection 



io|. 



To identify roles of 0o and it is convenient to consider the uniform case with no 
external magnetic field {V{x) = B = 0) and no supercurrent {q = 0). In this simple 
case, because of p{x) = po, Eq. (HM gives the gap less spectrum with quadratic dispersion 
Eq = p^/(2m). This means that 0o describes the ferromagnetic spin wave mode, which we 
call the transverse spin wave 1^. In addition, Eq. f llSp gives E_i = /(2m) +2|ci|po with a 
finite excitation gap 2|ci|po- This mode is called the quadrupolar mode [isl. The transverse 
spin wave described by 00 changes the total magnetization M by gpB- On the other hand, 
the quadrupolar mode described by 0_i changes M by 2gp-B. 

The current-flowing ferromagnetic state in the uniform case is stable only when the 
excitation energies of both the transverse spin wave 

Eo = ^-^ + 91^bB (16) 
2m 2m 

and the quadrupolar spin wave 

2 fc2 2 

^-1 = £^ - ^ + 2|ci|po + 2gpBB (17) 

are non-negative for arbitrary p. Since Eq < E^i, the stability condition of the current- 
flowing ferromagnetic state is given by 

^ < gf^BB. (18) 
2m 

In the following, we consider the cases where Eq. ( |T8|) holds. 



III. TUNNELING PROPERTIES OF TRANSVERSE SPIN WAVE 



To examine the stationary solution of the transverse spin wave described by (j)o{x,t), we 
set 4>o{x,t) = e~*'^*/Vo(a;) in Eq. (!T^ . We then have 



'2d^ 



+ v{x)-i-'^ + \^+i{x)\^ + n 



(j)o{x) 



(19) 



In obtaining Eq. f|T9|) . we have introduced the dimensionless variables, E = E/(c+po), 
X = x/i, q = g^, Ci = Ci/c+ (i = 1,2), V = V/{c+po), $+i = $+1/^, 0+i = 0+i/^, 
Q = qpbB/ (c+po)) and i = tc+po/h, and have omitted the bars from these rescaled variables, 
for simplicity. We always use these dimensionless units in the following discussions. 
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In the uniform system, Eq. (fT9|) has the plain wave solution (poi^) = e 
energy 

E 



='=''^^ with the mode 



1 



(20) 



Note that the stability condition Eq. ( ITSj) of the ground state is expressed as 2fi — > 
in the dimensionless units. For a given excitation energy E, the spin wave momentum far 
away from the barrier is given hyk= ^2{E + g2. 

We prove that perfect transmission of the transverse spin wave occurs at \k\ = g for a 
generic barrier potential V{x). With use of (|20i) . we first note that Eq. (1191) reduces to 

,2 







1 T./ X 



(j)o{x) 



(21) 



when k = q. When (polx) is replaced by the ferromagnetic condensate wave function 
$+i(x) = $+i(x, g) with the momentum q, Eq. coincides with the GP equation, 

,2 







-2^ + ^(^) 



l-|- + |$+i(x,g)p 



•^+1(3;, g), 



(22) 



which is nothing but Eq. ([8]) in terms of the dimensionless units introduced previously. We 
also note that when (po^x) is replaced by ^\i{x,q), Eq. ( 1211) coincides with the complex 
conjugation of Eq. f l^2]) . It immediately follows that $_|.i(x,g) and $^^^(2;,^) are linearly 
independent solutions to Eq. fl2T|) . The general solution to Eq. fl2Tl) is then expressed as 



0o(x) =a<l>+i(x, g) + /3$+i(a;, g). 



(23) 



where the coefficients a and /3 are determined so that 
of the tunneling problem 



satisfies the boundary conditions 



(x«-l), 
(a; > 1), 



(24) 



for k = q. When we set the phase of $+i(a;, g) so as to satisfy ^p+i{x = 0) = 0, the asymptotic 
form of the condensate wave function may be written as $+i(x, g) — )■ e'^[<i^+^s^i^)'^'fi/'^] ^ 
1) 26|, where Aip is the phase difference introduced in Eq. ffTOl) . From this asymptotic 
behavior of $+i(x, g) and the compatibility of Eqs. (123|1 and (12^ . one obtains (a,/3) = 
(e'^^/^O), i.e.. 



(/.o(x)=e*^'^/2^+i(x,g). 



(25) 
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and (Tq,Ro) = (e*^'^,0) at k = q. Perfect transmission |Top = 1 at = g for an arbitrary 
barrier potential V{x) thus follows. 

The same discussion is applicable to the case when k = —q. In this case, note that 
the condensate wave function in the supercurrent with the momentum —q (= $+1(2;, —q)) is 
related to the condensate wave function with the supercurrent momentum +q (= +q)) 
as —q) = +q). By repeating the above discussion by simply replacing $+1 by 

one may derive perfect transmission at A; = —q. 

The above argument applies when q = qc and hence the transverse spin wave transmits 
perfectly when q = qc and k = q^. This is contrasted with the transmission property of the 
Bogoliubov mode, which does not exhibit the perfect transmission when = Q'c js, 8, 10|. 

In the absence of supercurrent (g = 0), Eq. f l2^ does not necessarily imply that the 
wave function of the transverse spin wave mode {k = 0) has the the supercurrent behavior 
(i.e., the same form as the supercurrent wave function). The above derivation of perfect 
transmission does not apply because the incident and reflection waves are indistinguishable 
at = 0. When we discuss the tunneling problem in the limit A; — ?■ 0, we have to consider 
the problem with k being small but finite, and then take the limit /c — > 0. Following this 
procedure, we can show that this spin wave (with the momentum A; ^ 1) has the same 
form as the supercurrent wave function with the same momentum k with the accuracy of 
0{k) and from this fact, the perfect transmission follows in the limit k ^ when q = 0. 
The derivation is given in Appendix A. To conclude, the transverse spin wave, as well as 
the Bogoliubov mode, always has the supercurrent behavior, when the perfect transmission 
occurs. 

Let us see examples. Figure [3] shows transmission probability |Top of the transverse spin 
wave in the presence of the Gaussian- type barrier Eq. ([9]) with Vq = 2c+po- From the scaled 
plot (b), we see that the perfect transmission occurs aX k = q. 

Another example is provided by the exact solution in the presence of a 5-functional barrier 
V{x) = Vo6{x). In the case, the GP equation ([H]) has the following solution jsl, 

$+i(a;) = e'[''"-^§"(")^''][7(x) - sgn{x)iq], (26) 
where 7(0;) = a/1 — g^tanh[A/l — q'^{\x\ + xq)]- The phase factor 6q is given by 

^.e, ^ -Jl^ (27) 
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FIG. 3: Transmission probability |Top of the transverse spin wave, when the barrier is given by 
Eq. ([9]) with Vq = 2c+po- (a) momentum dependence, (b) |Top as a function of k/q. In obtaining 
this result, we have used the finite element method, imposing the boundary condition in Eq. (I24p . 
This method is also used to obtain Fig. [3 



where 7^ denotes \/\ — tanh(A/l — Q'^Xq) and Xq is determined from the boundary condi- 
tion on the derivative of at x = 0, which gives 



(1 - g^)7g - 7g 
72 + g2 



ikx 



(28) 



(29) 



Substituting Eq. ( I26l) into Eq. ( lT9l) . one obtains 

0o(x) = [— sgn(x)iA; + 7(x)]e' 
To solve the tunneling problem, we set 

0o(x) = + 7(x)]e*'''' + RQ\-%k + 7(x)]e"*''^, (x < 0), 
= To[-zfc + 7(x)]e+^'=^, (x > 0). 

As usual, the coefficients Rq and Tq are determined so as to satisfy the boundary conditions. 



(30) 



0o(+O) = 0o(-O), 



c?0o(x) 
dx 



x=+0 



#o(x) 
dx 



(31) 



z=-0 



11 




FIG. 4: (a) Transmission probability \Tq\^ of the transverse spin wave in the case of a (5-functional 
barrier. Panel (b) also shows |Top as a function oik/qto clearly show that the perfect transmission 
always occurs &i k = q. We take Vq = 2{cj^n)/^. 



1 



We then have 

To =ik 

Ro = -l + ik 



7g 



+ -q^ + k-fg{i + k-fg) + q^{l + k^ --f^) 

1 , % 



-ik + 7q -q^ + k^q{i + k^q) + q'^{l + A;^ - 7^ 



(32) 
(33) 



In the absence of a superflow (g = 0), Eq. fl32|) reduces to 

^7o(l + k') 



Tn 



(34) 



Thus, perfect transmission (|Top = 1) is obtained in the long wavelength limit (A; — )■ 0) 12 



14 



16(1 . In contrast, Eq. ( 132|) vanishes in the limit A; — when g is finite. This shows clearly 



that, in contrast to the Bogoliubov mode, anomalous tunneling of the transverse spin is not 
realized in the long wavelength limit in the presence of a finite superflow. 
When /c = g > 0, Eq. fl32|) becomes 



Tn 



(35) 



which gives perfect transmission |Top = 1. In the scalar BEC, it has been pointed out 
that the wave function of the Bogoliubov mode has the same form as the condensate wave 
function when the anomalous tunneling occurs in the absence of supercurrent jsl and in the 
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presence of supercurrent 8| . In the present case, from the comparison of Eq. ( 126|) with Eq. 
fl30|) at k = q, one finds 

0o(a:) = e-*^^$+i(x,g), (36) 



as in the case of generic potential barrier fl25|) . The fact that perfect transmission of the 
transverse spin wave occurs at the finite momentum k = q is different from the condition of 
anomalous tunneling in the case of the Bogoliubov mode; however, the supercurrent behavior 
can be also seen in the transverse spin wave when perfect transmission occurs {k = q). In 
the absence of supercurrent (g = 0), we show in Appendix B that this spin wave (with the 
momentum k <^ 1) has the same form as the supercurrent wave function with the same 
momentum k with the accuracy of 0{k). This fact illustrates the supercurrent behavior of 
the transverse spin wave in the limit A; — ?■ when q = 0. 

We briefly note that the amphtude transmission coefficient in Eq. ( |32l) satisfies To(— fc) = 
To{k)*. Because of this symmetry, we find that perfect transmission also occurs at k = —q. 

Figure m shows the transmission probability |Top of the transverse spin wave through the 
5-functional barrier. When the supercurrent momentum q becomes finite, panel (a) shows 
that the value of k at which perfect transmission (|Top = 1) is realized becomes large. At 
perfect transmission, panel (b) confirms that \k\/q = 1 is satisfied, as expected. Figure H] 
also shows that the perfect transmission occurs only when \k\ = q. 



IV. TRANSMISSION PROPERTIES OF QUADRUPOLAR MODE 



The quadrupolar spin wave is described by the fiuctuation component which obeys 
Eq. (fT5|) . In the stationary state, 0_i(x, t) = 0_i(x)e~*^*, this equation reduces to, in the 
dimensionless form, 



1 d 



q^ 



2^ + nx)-l-^ + (l-2c,)|$.. 



2n 



(37) 



The mode energy in the uniform system is obtained as E = {k"^ — q'^)/2 + 2|ci| + 2Q, which 
is positive for arbitrary k when the stability condition (1181) {2Q > q^ in dimensionless units) 
of the current-fiowing ground state is satisfied. Under this condition, the momentum k has 
the form k = ^/2[E + 2(ci — Q)] + g^. In considering the tunneling problem, the boundary 
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FIG. 5: Transmission probability jTLip of the quadrupolar mode. The barrier potential is given 
by Eq. Q with Vq = 2c+po. The ratios ci/c+ of the coupling constants are (a) -0.01, (b) -0.1, and 
(c) -1. 

conditions far away from the barrier are given by 

4>-i{x) = exp (ikx) + exp {—ikx), {x <^ —1), 

= r_i exp (tkx), {x > +1). (38) 

Figure shows the transmission probabiUty |T_ip of the quadrupolar mode, when the barrier 
potential is given by Eq. In contrast to the Bogoliubov mode and the transverse spin 
wave (which are described by Eqs. (fT3|) and (JH]), respectively), the interaction effects on 
the quadrupolar mode cannot be simply described by the single parameter c+ = Ci + C2. 

14 



Because of this, the couphng ci remains in the dimensionless equation ( 137|) . (Note that ci 
in Eq. fl37|) actually denotes ci = ci/c+.) Indeed, we find in Fig. [5] that the transmission 
probability strongly depends on the value of Ci. Although the quadrupolar spin mode could 
show perfect transmission at a certain value of k (See, for example. Fig. [S] (a)), it owes 
to a fine tuning of the parameter. The tunneling property of this mode is rather close to 
the ordinary quantum mechanical tunneling of a particle in the sense that perfect refiection 



occurs when k 
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V. EFFECTS OF SPIN-DEPENDENT BARRIER 

So far, we have considered the case when a barrier potential couples to the local density of 
particles as V{x)p{x). Since the transverse spin wave is a Nambu-Goldstone mode associated 
with spontaneous breakdown of rotational spin symmetry, the tunneling property of this 
mode is expected to be strongly affected by a spin-dependent magnetic barrier. To see this, 
we examine tunneling through a magnetic potential barrier 

V,{x) = -J{x)F,{x). (39) 

Although this tunneling problem is somehow academic in cold atom gases, it is still in- 
structive to see how the perfect transmission of the transverse spin wave is related to the 
symmetries broken or preserved by the barrier. For simplicity, we assume the absence of 
supercurrent (g = 0), as well as the external magnetic field (5 = 0). We also take J(x) > 0, 
which enhances the Sz = +1 spin component near the barrier. When we replace V{x)p{x) 
in Eq. by Vs{x), we again obtain a 3 x 3-matrix GP equation similar to Eq. (E]). By nu- 
merically solving the GP equation with the boundary condition ($+i, $0) '^'-i) = (v^' 0) 0) 
for |x| ^ 1, we have confirmed that $o = '^'-i = even near the magnetic potential Vs{x). 
Thus, one may consider the single-component GP equation for $_|_i(x), given by 

1 







2d.- + 



(40) 



Equations for the collective modes can be also obtained in the same manner as we have 
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done in Sec. II. The results are 



'+1 



dt 



dt 



dt 



- J{x) - 1 + 2p{x) 



'+1 



1 . 



+ J(s)-l + (l-2ci)p(x) 



(41) 
(42) 
(43) 



Equations (HOj) and (14T!) have the same forms as the GP equation and the Bogohubov 
equation in the scalar BEG, respectively. Thus, the Bogoliubov mode exhibits perfect trans- 
mission through the magnetic barrier Vs{x) in the low energy limit. 

To see whether or not the transverse spin wave shows perfect transmission, it is convenient 
to consider the 5-functional barrier J{x) = Jo6{x) (Jo > 0). In this simple case, the GP 
equation (HOl) has the exact solution 



$+i(x) = coth (|x| + Xq 



where Xq = coth ^[(a/Jq + 4 — Jo)/2]. Substituting p{x) 
finds that the transverse spin wave solution has the form 



(44) 

|$+i(x)p into Eq. (02]), one 



(f)o{x) = f^{k,x) + Rof^{-k,x), (x<0), 
(f)oix) =Tof+{k,x), ix>0), 



(45) 



where f±{k,x) = [^ik + coth(±x + Xo)]e*'^^. Since the potential barrier J{x) is absent in 
Eq. ( l42l) . the boundary conditions at x = are simply given by the matching of the wave 
function and its gradient, which give 

fc(l + fc2) 



To 

Rq 



{k + iT){-l + k^ + ikT + T^)' 
-iPfl -r2) 



(46) 
(47) 



where F = coth(Xo 



{k + tT){-l + k^ + ikT + T^)' 

The coefficient Tq in Eq. P7|) vanishes in the limit — )■ 0, unless F = 1 
(i.e., Jo = 0). Thus, in contrast to the 'non- magnetic' potential V{x) discussed in Sec. Ill, 
the magnetic barrier Vs{x) leads to perfect reflection of the transverse spin wave in the low 
energy limit. 

Figure |6] shows the transmission probabilities of the three collective excitations through 
the magnetic barrier potential. As discussed above, the Bogoliubov mode shown in panel 
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FIG. 6: Transmission probability through the magnetic potential barrier when J{x) = Joe~^ . 
(a) Bogoliubov mode, (b) transverse spin wave, (c) quadrupolar spin wave. The inset shows the 
spatial variation of the ferromagnetic condensate wave function <I>+i(x). 

(a) only exhibits the perfect transmission at /c = 0. Although the transverse spin wave does 
not feel the magnetic potential (See Eq. (H2l) ). it is still scattered by the 'barrier' formed by 
the condensate wave function around x = shown in the inset of panel (a). Since Eq. ( H2ll 
has the same form as the ordinary one-particle Schrodinger equation with the 'potential 
barrier' — 1, the perfect reflection at k = can be understood as a usual property 

in one-particle quantum mechanics. 

As we have seen, the spin-dependent barrier suppresses the perfect transmission of the 
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FIG. 7: Transmission probability of the Bogoliubov mode through the barrier given by Eq. (j48p . 
We take G{x) = Goe-^^/^\ 

transverse spin wave but the Bogoliubov mode still exhibits anomalous tunneling. This 
fact implies that the perfect transmission (or the anomalous tunneling phenomenon) of a 
Nambu-Goldstone mode is strongly affected by the barrier that breaks the symmetry related 
to the origin of the mode. To see this in a simple manner, we replace V{x)p{x) in Eg. ([3l) 
by the one that couples to the local phase (p+i{x) of the condensate wave function j27|, 
such that 

Vg{x) = G{x) \<i>_^^(^x)\ " cos(^+i(x). (48) 

When G{x) = GoS{x), the transmission probability jT+ip of the Bogoliubov mode with the 
incident momentum k is given by 

(We summarize the derivation in Appendix C.) Here, r] = A;[|Go|^^ + '^(^^ + '*^)]) and 
K = ^2{y/l + 1), where E = {k / y/2) y^{kpI^^^T2 is the mode energy. Equation (|49|1 
vanishes in the low energy limit — ?■ 0, so that perfect reflection occurs. Even for a barrier 
with a finite width, this conclusion is unaltered, as shown in Fig. [71 (Equations to obtain 
this figure are summarized in Appendix C.) The above two results for the transverse spin 
wave and the Bogoliubov mode indicate that the anomalous tunneling phenomenon is deeply 
related to the invariance of the barrier potential with respect to continuous symmetries that 
are broken spontaneously, such as the U(l) gauge and spin rotation. 
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VI. SUMMARIES AND DISCUSSIONS 



To summarize, we have investigated tunneling properties of the collective Bogoliubov 
mode, as well as spin wave excitations (transverse spin wave and quadrupolar mode), in 
the ferromagnetic state of a spin-1 spinor Bose-Einstein condensate. Within the framework 
of the mean-field theory at T = 0, we have calculated the tunneling probability across a 
barrier potential of each mode as a function of the incident momentum, to see whether or 
not anomalous tunneling phenomenon occurs. 

We showed that, as in the case of a scalar EEC, the Bogoliubov mode also exhibits the 
anomalous tunneling phenomenon in the ferromagnetic state of a spin-1 spinor EEC in the 
presence of a barrier that couples only to local density. Namely, perfect transmission always 
occurs in the low energy limit, even in the presence of supercurrent, unless the superfiow 
is in the critical supercurrent state. In the critical supercurrent state, tunneling of the 
Bogoliubov mode is accompanied by a finite reflection. All these results are the same as the 
case of a scalar EEC 

The transverse spin wave, which is the other Nambu-Goldstone mode associated with the 
broken spin rotational symmetry, also shows the anomalous tunneling phenomenon in the 
presence of a barrier that couples only to local density. However, while this phenomenon 
occurs in the long wavelength limit in the case of the Bogoliubov mode, perfect transmission 
of the transverse spin wave is realized when the momentum of the spin wave coincides with 
that of supercurrent. In addition, this phenomenon still occurs in the critical supercurrent 
state, which is also to be contrasted with the case of the Bogoliubov mode. On the other 
hand, when perfect transmission occurs, the wave function of the transverse spin wave has 
the same form as the condensate wave function in the current carrying state. This property 
is shared with the Bogoliubov mode in the long wavelength limit. 

This result is important in the context described below. The perfect transmission of the 
transverse spin wave in the current-flowing state provides much simpler example of anoma- 



1^ 



lous tunneling, compared to earlier cases. In the earlier cases (the Bogoliubov mode 

nn nnn 

I 111 . |13| , the transverse spin wave in the absence of supercurrent [12|, |14|-|16| ) , the anomalous 
tunneling always occurs in the limit where the wavenumber k approaches zero. In the present 
case, on the other hand, anomalous tunneling occurs at a finite wavenumber. When we dis- 
cuss the tunneling problem in the limit A; — )■ 0, we have to calculate the transmission rate 
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for finite k first, and then consider the limit A; — )■ 0. This limiting procedure is unnecessary 
when we consider the tunneling problem of particles with finite wavenumber. This is the 
reason why the proof of perfect transmission for g 7^ presented in Sec. II is simpler than 
those for g = given in [15] and Appendix A in the present paper. Owing to this simplicity 
of the perfect transmission of the spin wave in the current-flowing state, we can prove the 
perfect transmission of the spin wave for a general potential barrier and g 7^ as a direct 
consequence of the coincidence of wave functions between the spin wave and the flowing 
condensate. 

We have also examined tunneling properties of the quadrupolar mode. This collective 
mode has a finite excitation gap, and is not the Nambu-Goldstone mode. In contrast to the 
above two Nambu-Goldstone modes, the quadrupolar mode does not exhibit the anomalous 
tunneling behavior. 

The tunneling properties of excitations are determined by the properties of potentials 
as well as characters of excitations. When a barrier potential locally destroys a symmetry 
that is spontaneously broken in the spinor BEG, the anomalous tunneling phenomenon 
of the corresponding Nambu-Goldstone mode disappears. Indeed, perfect transmission of 
the transverse spin wave does not occur, when the barrier potential depends on spin and 
locally breaks spin rotation symmetry. Similarly, perfect transmission of the Bogoliubov 
mode does not occur either, when the barrier potential depends on the local phase of the 
condensate wave function and locally breaks U(l) gauge symmetry. The presence of the 
symmetry breaking potential inhibits the wave function of a Nambu-Goldstone mode from 
coinciding with the condensate wave function. These results indicate that the symmetry 
that the barrier potential possesses is crucial for the anomalous tunneling phenomena of 
Nambu-Goldstone modes. We infer that anomalous tunneling occurs in a Nambu-Goldstone 
mode only in the presence of symmetry-preserving barrier potentials. This result shown in 
Sec. |V]is significant because this is the first study on relation between anomalous tunneling 
and symmetry-property of potentials. 

The present paper has treated the tunneling problem for all the three excitations sep- 
arately. It seems to be difficult to create an excitation as a single mode from the spin-1 



BEG to examine tunneling properties of excitations experimentally; but, this difficu 



be overcome by making use of the Bragg scattering and a strong magnetic field 



ty can 



lej. The 



Bragg scattering can selectively create an excitation with a certain momentum p and a cer- 
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tain energy E. When energies of the three collective excitations at small momenta are well 
separated from each other by using the Zeeman effect under a strong magnetic field, one 
would be able to selectively create one of the three excitations as a single mode. 

Since low energy properties of the ferromagnetic spinor BEC are dominated by the three 
collective excitations, our results would be useful for understanding dynamical properties in 
this superfiuid phase. 

In this paper, we have focused on the ferromagnetic phase of a spin-1 BEC. It is an 
interesting problem to extend our present work to the polar phase, which is the other phase 
of the spin-1 BEC. We have examined the same problem for the spin wave modes in the 



current carrying state of the polar phase 28|. We found that the anomalous tunneling 



phenomenon of the spin wave modes occurs unless the system is in the critical current state. 
The anomalous tunneling phenomenon in itself is the same character as the Bogoliubov 
mode and the transverse spin mode in the ferromagnetic phase. However, in the critical 
current state, the tunneling property is qualitatively changed by the interaction parameter. 



These findings will be reported in a separate paper 
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Appendix A: Proof of perfect transmission of transverse spin wave for a generic 
barrier at q = 

The proof of perfect transmission of the transverse spin wave for g 7^ is presented in 
Sec. mil but this proof is not applicable when q = 0. While the proof for the case with g = 
has been given in 15|, we present in this appendix another proof of perfect transmission of 
the transverse spin wave at g = in order to show that the similarity between the wave 
function of an excitation and the supercurrent wave function leads to perfect transmission 
for the case with g = as well as the case with g 7^ 0. In this appendix, we abbreviate 
g) as g) for simplicity. 
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First we introduce 

f,{x) = mx,q)\^-mx,0)\^ (Al) 

which is Oi^q"^) because depends on q only through (Equation for g)p 

contains q only through ). When q = 0, the equation of motion for the transverse spin 
wave with the momentum A; > 



is rewritten as 



with 



hk - fk{x) 00 (a;) = 



(x) = 



1 



(A2) 

(A3) 
(A4) 



hk ^ + V{x) - 1 - y + k)\'. 

We seek for the solution to ( 1A2I) satisfying the boundary condition Eq. (IMl) for the tunneling 
problem, assuming the form 



with a normalization factor A^^. Equation (1A3P is rewritten as 



(A5) 



hkA(l)k{x) = ^fkix)(j)o{x). 



(A6) 



Substituting the following expressions: 



AMx) = Xk{xMx,k) (AT) 

PX POO 

Xk{x) = 2N^^ dx'[^x\k)Y^ dx"fkix")Mx"Mx",k) (A8) 

Jo Jx' 

into flA6p . we see that flA7p with flASp yields a solution that we seek for. However, it is 
a formal solution; Equation (lASp contains 0o in the right-hand side. In the following, we 
consider the asymptotic form oi Xk{x) at \x\ 1. We first note that |$(x — > ±oo, A;)p — )■ 1 
irrespective of /c and hence fk is localized around a; = 0. We then see that 

(A9) 



dx"fk{x")Mx"Mx",k) 
approaches zero when x' — )■ oo while it approaches a constant 



At = 



dx"fkix")M^"Mx",k) 



(AlO) 
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when x' — oo. From this property and the asymptotic behavior $(x — )■ ±00, A;) — )■ 
gj[A:x+A(/PfcSgn(x)/2] q£ ^^iq condensate wave function, it follows that 



(x > 11 



-Afc(e-2^^^ - l)e'^^>'/{ik) + 6^, {x < -1). 



(All) 



Here and in the following, we make /c-dependence of the phase shift A(p explicit. The 
constants a^, bk are, respectively, given by 



/•oo POO 

ak = 2 dx'i^kix')]-^ / dx"fi 

Jo Jx' 



"Mx")Mx")<!>{x",k), 



bk = -2 I dx' 

-00 







J ~oo J —00 

Substituting (IX7|) with flMj) into flXSj) and using flAlip . we obtain 



dx"fk{x")Mx")<^{x",k) 



(A12) 



6o(x) ^ 



[Nk - i{Xk/k)e'^'^^ + bk] e-'^^^/^e''''' 



(A13) 



From flAlSp and ([21]), the solution with 



iV, = e^W-6fc + z(Afe/A;)e*^^S 



To = e 



+ z{Xk/k)e'^^^''/^ + {ak-bk)e'^'^'/\ 



Ro = i{Xk/k)e 



(A14) 
(A15) 
(A16) 



follows. 

As noted previously, fk = 0{k'^) when \k\ ^ 1. (I)q{x) is at most 0{1) (when the perfect 
reflection occurs in the limit /c — ?■ 0, the wave function is 0{k)) and $(x, /c) is 0{1). The 
quantities a^, bk-, Xk contain commonly fk(f)o^{x, k) as a factor in their integrands and hence 
they are (9(/c^). The phase shift Aipk is 0{k) for \k\ <^ 1. We thus obtain 



To = 1 + zA; lim (A(^fc/fc + Xk/k^) + 0{k^), 



Ro = tk lim(Afc/A;2) + C(A;2) 



(A17) 
(A18) 



and the perfect transmission |Top — )■ 1 occurs in the limit — 0. 
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Appendix B: transverse spin wave solution in the low energy limit in the absence 
of supercurrent 



In this appendix, we show that the transverse spin wave mode has the same form as 
the current carrying condensate wave function in the long wavelength limit, even when the 
supercurrent is absent (g = 0) and the potential barrier is given by V{x) = Vo6{x). 

Assuming that the supercurrent momentum q is small, we expand the condensate wave 
function in Eq. fl26|) in terms of q up to 0{q), which gives 

,i[gx-sgn(x)e,]^-^^^ - zgsgn(x)] = 



(Bl) 



Here, 7(x) = tanh(|x| + Xq), and 9g = —q/'y{0). We also expand the spin wave solution in 
Eq. f l30p around k = when q = 0. The amplitude transmission and reflection coefficients 
in Eqs. and are expanded as, respectively, 

^ ' (B2) 

(B3) 



Ro =ik 7(0) 



7(0)7 • 



We then have, within the accuracy of 0{k), 

0o(x<O) = [+zA; + 7(x)]e*'=" + ^A; (^7(0) 

[7(x) + ik] 



ikx 



7(0) 
7(0;) + ik 



ik + 7(x)]e' 



-ikx 



(B4) 



and 



> 0) =Hfc + 7(x)]e*'^^' 



1 + 7(0) + - 



7(0) 



(B5) 



Comparing Eqs. flB4p and fIBSP with Eq. fIBll) . we have 



Jk^{0) i[kx-~sgn{x)ek] 



[y{x) — iksgn{x)] + 0{k'^ 



=e*'=^(°)$+i(a;,A;) + C(/t2). 



(B6) 



Thus, the wave function of the transverse spin wave in the long wavelength limit [k — )■ 0) has 
the same form as the condensate wave function in the supercurrent state with the momentum 
k within the accuracy of 0{k). 
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Appendix C: Tunneling through a barrier coupUng to the phase of condensate 



In the case of Eq. (HHj) . the GP equation is given by 

1 + p{x) + G{x) 



2dx^ 



^+i(x) 



0. 



(CI) 



2|$+i(x)|3 

In this case, when G{x) < {G{x) > 0), the uniform solution = 1 ($+i = — 1) gives 

the minimum energy. When G{x) = Gq6{x), the Bogoliubov equation for the Bogohubov 
mode is then given by 

M<l>-,iix,t) ^ ^ ^^^^^^^^ _ ^^^^^^^^ ^ ^\Go\S{x)[M^,t) + <Pl,ix,t)]. (C2) 

As usual, the last term in Eq. ( ]C2I) may be conveniently absorbed into the boundary condi- 
tions at X = 0. Setting 0+i(x, t) = u{x)e~^^^ — f *(x)e*^*, one obtains an eigenvalue problem 




~2dx^ 



+ 1 -1 

1 d'' 



\ 



2 c/x2 V 




The basis functions of general solutions to ( 1C3|) are given by 





Here we introduce k = \ 2{yjl + - I), k = \ 2{Vl + + 1), and 



A 



\ 



1 / VTTe^ 



+ 1 



B 



\ 



1 / VTTe^ 

2 \ E 



which are normalized so as to satisfy — B"^ = 1. 

Under the condition that the Bogoliubov mode is injected from x = — oo, we set 

(x < 0), 
(x > 0). 

The boundary conditions at x = are the matching of Eqs. (1060 and (ICTO . and 

I Go I 


















1 e-^^ 






x=+0 




x=-0 



u{0)+v{0) 
u{0)+v{0) 



(C3) 



(C4) 



(C5) 



(C6) 
(C7) 



(C8) 
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The coefficients in Eqs. (1C6|1 and (1C7P are determined from these boundary conditions as 

V 



(Cll, 

where ri = kllGolk"^ + nik"^ + /t^)]. The transmission probabihty (Eq. ( I49l) ) follows from 
Eq. (|C9|. 

When the barrier has a finite width (e.g., G'(x) = Goe~^^^^^), the Bogoliubov equation, 
corresponding to Eq. flC3l) . becomes 




1 + 2p(x) + xix) 



( iil 



1 rf2 

-2^ - 1 + 2^>(-) 




where 



Xix) =^G[x)- 
W{x) =\g{x\ 



1 



1 



<l>*(x) 



2<l>fx) 



3 |$(x)p- 
2 <l>(a; 
3|$rx)P-$2fx' 



(C12) 

(C13) 
(C14) 



2 <l>*(x) 

In this case, we numerically solve Eq. (IC12I1 . regarding Eqs. (IC6I1 and ( 1C7I1 as the boundary 
conditions at x = — oo and x = oo, respectively. 

Before closing this appendix, we briefly remark on transmission of the transverse spin 
mode through the potential ( HHj) . which couples to the local phase of the condensate wave 
function and gives the uniform solution of In this case, the equation of the 

excitation for 5"^ = is given by ([H]), where p{x) is spatially constant (i.e., p{x) = po) 
and V{x) = 0. (The potential ( l48l) does not enter in ( fT4l) . since it does not couple to the 
component Sz = 0). As a result, the transverse spin mode propagates through a uniform 
medium, and shows perfect transmission. 
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